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Abstract 



Close binary systems of compact objects with less than ten minutes remaining 
' before coalescence are readily identifiable sources of gravitational radiation 

for the United States Laser Interferometer Gravitational-wave Observatory 



(LIGO) and the French-Italian VIRGO gravitational-wave observatory. As 



a start toward assessing the full capabilities of the LIGO/VIRGO detector 



network, we investigate the sensitivity of individual LIGO/VIRGO-like inter- 
ferometers and the precision with which they can determine the characteristics 
of an inspiralling binary system. Since the two interferometers of the LIGO 
bJQ' detector share nearly the same orientation, their joint sensitivity is similar to 

that of a single, more sensitive interferometer. We express our results for a 
single interferometer of both initial and advanced LIGO design, and also for 
' the LIGO detector in the limit that its two interferometers share exactly the 

same orientation. 

We approximate the secular evolution of a binary system as driven exclu- 
sively by its leading order quadrupole gravitational radiation. Observations of 
a binary in a single interferometer are described by four characteristic quan- 
tities: an amplitude A, a chirp mass A4, a time T, and a phase ip. We find 
the amplitude signal-to-noise ratio (SNR) p of an observed binary system as 
a function of A and A4 for a particular orientation of the binary with respect 
to the interferometer, and also the distribution of SNRs for randomly ori- 
ented binaries at a constant distance To assess the interferometer sensitivity, 
we calculate the rate at which sources are expected to be observed and the 
range to which they are observable. Assuming a conservative rate density 
for coalescing neutron star binary systems of 8 x 10~ 8 yr™ 1 Mpc~ 3 , we find 
that the advanced LIGO detector will observe approximately 69 yr -1 with an 
amplitude SNR greater than 8. Of these, approximately 7 yr" 1 will be from 
binaries at distances greater than 950 Mpc. 
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We give analytic and numerical results for the precision with which each 
of the characteristic quantities can be determined by interferometer observa- 
tions. For neutron star binaries, the fractional one-sigma statistical error in 
the determination of A is equal to l/p. For p > 8, the fractional one-sigma 
error in the measurement of M in the advanced LIGO detectors is less than 
2 x 10~ 5 , a phenomenal precision. The characteristic time is related to the 
moment when coalescence occurs, and can be measured in the advanced de- 
tectors with a one-sigma uncertainty of less than 3 x 10 -4 s (assuming p > 8). 

We also explore the sensitivity of these results to a tunable parameter in 
the interferometer design (the recycling frequency). The optimum choice of 
the parameter is dependent on the goal of the observations, e.g., maximiz- 
ing the rate of detections or maximizing the precision of measurement. We 
determine the optimum parameter values for these two cases. 

The calculations leading to the SNR and the precision of measurement 
assume that the interferometer observations extend over only the last several 
minutes of binary inspiral, during which time the orbital frequency increases 
from approximately 5 Hz to 500 Hz. We examine the sensitivity of our re- 
sults to the elapsed time of the observation and show that observations of 
longer duration lead to very little improvement in the SNR or the precision 
of measurement. 

PACS numbers: 04.80.+z,04.30+x,97.60.Jd,06.20.Dk 
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I. INTRODUCTION AND MOTIVATION 



Both the United States Laser Interferometer Gravitational-wave Observatory (LIGO 
0J3D and the French/Italian VIRGO gravitational- wave observatory || are expected to 
begin operation in the late 1990s. Inspiralling binary systems of compact objects — either 
neutron stars or stellar mass black holes with orbital frequencies ranging from 5 s^ 1 to 
500 s _1 — are currently regarded as the most certain observable source for these detectors: 
the density of sources [^]|| suggest event rates of several per year and the radiation they 
emit can be calculated unambiguously j6|,|7|j8|j9|,p!0|, p!T| , p^] . In support of the LIGO/ VIRGO 



observational effort, theoreticians must combine refined calculations of the radiation from 
binary inspiral with the anticipated detector properties to deduce the instrument sensitivity. 
These calculations will, in turn, play an important role in the final design and ultimate use 
of the instruments. In this paper, we begin a detailed analysis of binary inspiral as a source 
of gravitational radiation for a LIGO/VIRGO-like interferometric detector. 

Our goal in this work is to estimate the sensitivity of LIGO / VIRGO-like interferometers 
and the LIGO detector by determining the rate at which inspiralling binary systems can be 
detected, the range to which they can be observed, and the precision with which they can 
by characterized in a single interferometer. We also explore the compromises that must be 
made as different questions are asked of the observations. The optimal design and operation 
of these interferometers depends on a detailed understanding of the nature of the detector 
response to the radiation, the detector noise power spectral density (PSD), and the questions 
the observation is meant to resolve. For example, we show in §[V] that the goal of observing as 
many sources as possible (without necessarily being able to characterize them precisely) leads 
to a different optimal interferometer configuration than the goal of characterizing observed 
sources as precisely as possible (while allowing that weak sources may be missed entirely). 

The ultimate goal of our assessment of binary inspiral is to determine the ability of a 
single interferometer (or a network of such interferometers) to 

1. Detect the gravitational radiation from the last few minutes of inspiral of a binary 
neutron star or black hole system, and 

2. Measure the parameters describing the detected binary system from the observed 
gravitational radiation. 

By detection we mean the determination of the presence or absence of a signal characteristic 
of an inspiralling binary system in the output of a detector, irrespective of the particular 
parameters that might characterize the observed binary system. By measurement we mean 
the determination of the parameters that characterize the signal presumed to be present 
in the detector output. In a real detector noise can mask or distort a signal present in 
the detector output; alternatively, it can conspire to appear as a signal characteristic of a 
binary system. Consequently, any conclusion we draw from observations (e.g., that we have 
detected an inspiralling binary system) is associated with a probability that characterizes 
our certainty in its validity. 

As a practical matter, reliable detection of gravitational radiation will initially require 
coincident observation between two or more interferometers so that non-Gaussian noise 
events can be differentiated from gravitational radiation signals. When completed, the 
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LIGO detector will consist of two interferometers: one in Hanford, Washington and one in 
Livingston, Louisiana || . The relative orientation of the two interferometers has been chosen 
to maximize their sensitivity to a single polarization state of the gravitational radiation 
impinging on the Earth! ||13|| ; consequently, the two interferometers act together like a 
single interferometer whose sensitivity will be greater than that of either of its components. 
This analogy is not exact: combining the interferometers in this way ignores the different 
arrival times of the signal at the two distinct interferometers. Nevertheless, our single- 
detector analysis is directly relevant to the actual LIGO configuration, and we discuss our 
results in the context of (i) single LIGO/VIRGO-like interferometers of the proposed initial 
and advanced designs, and (ii) the LIGO detector in the limit that its two independent 
interferometers share exactly the same orientation and ignoring the difference in the signal 
arrival times. 

In principle, simultaneous observation of an inspiralling binary system in three interfer- 
ometers of different orientations is sufficient to measure (among other characteristics) the 
binary's luminosity distance di, its position in the sky, and a function M. that depends 



only on the masses of its components and its cosmological redshift 0|T^] • If the position is 
known through other observations, then A4 and di can be measured from observations with 
only two interferometers of different orientations. Observation of an inspiralling binary in a 
single interferometer can measure M. and a waveform "amplitude" A, which depends on the 
luminosity distance and orientation of the binary with respect to the detector. Even this 
limited information is of astrophysical significance, however: from observations of the distri- 
bution of A and M. among the observed binaries one can determine the distribution of the 
component masses of inspiralling binary systems, find the Hubble constant and deceleration 



parameter, and test cosmological models [16 



The study of binary systems as sources of gravitational radiation began in 1963 when 
Peters and Mathews || made the first detailed calculation of the gravitational radiation 
luminosity from inspiralling binary systems, focusing on the leading order quadrupole ra- 
diation from two point particles in circular and elliptical orbits. Clark and Eardley [17] 



explored how gravitational radiation (among other effects) drives the orbital evolution of 



binary systems. Clark [18] suggested inspiralling binary systems as an important source of 



gravitational radiation for modern interferometric detectors; however, he considered only the 
burst of radiation from the coalescence event itself. Thorne recognized the importance of the 
gravitational radiation from the final few minutes of inspiral. Since 1987 a number of inves- 
tigators have worked with increasing sophistication on problems related to the observation 
of inspiralling binaries in interferometric detectors. Some have focused on understanding 
and refining interferometer detector technology [ 19 , 20| , pT| , p2| , p3] , pi]] , others have focused on 



the data analysis problems of detecting or determining the characteristics of a binary system 
from the radiation [^,^,^7j,^,^,Q , and still others have focused on refining our under- 
standing of the gravitational radiation waveform from these systems |7J,|3J0,ro,O,O,Rl|,R3,R3 



1 While primarily sensitive to a single polarization state of the radiation field, LIGO will, with 
more limited sensitivity, be capable of observing both polarizations simultaneously. In a subsequent 
paper we will take the exact detector orientations into account in our analysis. 
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or their rate of occurrence in the Universe ||||. 

The remainder of this paper is organized as follows. In §|| we review how the precision 
with which the parameterization of a signal, observed in a noisy detector, is determined. In 



[II| we apply these techniques to the particular problem of finding the precision with which 



the parameters of an inspiralling binary system can be measured. These results depend only 
on the detector noise power spectral density (PSD) and not on the type of detector (e.g., 
interferometer or bar). In § P"V| we apply the results of §[ITJ to both the initial and advanced 
proposed LIGO/VIRGO-like interferometers. We discuss the astrophysical implications of 
these results in §M and present our conclusions in Wl . 



II. MEASUREMENT AND UNCERTAINTY 



In this section we review the techniques we use to determine the statistical uncertain- 
ties in observations of binary inspiral. A complete discussion is found in Finn []34}| . The 
techniques developed there are closely related to those associated with signal analysis by 
optimal filtering. For more information on optimal filtering and signal analysis we direct 
the reader to Oppenheim, Willsky, and Young , the review by Davis |!6| and references 
therein, Wainstein and Zubakov ||37||, and Hancock and Wintz [j38fl. 



A. Introduction 

Consider a deterministic (i.e, not stochastic) source of gravitational radiation (e.g., an 
inspiralling binary system) and a detector (e.g., a laser interferometer). We write the re- 
sponse of the detector to the radiation as a superposition of noise n(t) and signal m(t; /it), 
where n is a minimal set of parameters that uniquely characterizes the detector response 
(absent the noise) to the radiation for the entire duration of the observation (in §JTJ we show 
that for an inspiralling binary system observed by a single interferometric detector /x is a 
four dimensional vector). Let the source of radiation be characterized by fi. In analyzing 
the output of the detector, we have two goals: 

1. To determine whether a signal is present in the detector output, and 

2. To find the precision with which we can determine jx (assuming a signal is present). 

Owing to the detector noise, we cannot determine with certainty either the presence 
of a signal or (assuming it present) ft. Instead, we find the probability that the detector 
output is consistent with the presence of a signal, and represent our uncertainty in ft by a 
set of volumes V(P) in parameter space, such that ft is in V(P) with probability P. The 
volumes V(P) are a measure of the sensitivity of the detector. Throughout this section, 
we will assume that we have determined that a signal is present so that the probability 
P associated with V(P) is a conditional probability. In a later paper, we will discuss the 
determination of the probability that a signal of the form m(t; fi) is present in the detector 
output. 

Each observed signal is immersed in its own realization of the detector noise; consequently 
we cannot know in advance of an observation what the volumes V(P) will look like. We can, 
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however, determine what the most likely volumes are for a given observation, and that is 
what we do here. In this section we review briefly the procedures used to find the most likely 
volumes V(P). These procedures are developed in Finn P^| , and we refer the interested 
reader there for more information. 



B. The probability volumes V(P) 



We characterize the observed output of our detector as a time series g(t), which is a 
superposition of noise n(t) and (perhaps) a signal m(t; p): 



g(t) 



n(t) + m(t; fx) signal present 
n(t) signal absent. 



(2.1) 



The parameter fx is fixed but unknown, and we assume the noise is drawn from a stationary 
process. The probability density that a signal with parameterization 11 is present in the 
detector output g(t) is 



where 



(The conditional probability that a \ 
signal characterized by \x is present | 
given the detector output g(t). 

Hp) 



P(0) = 

P{m) = 
A = 

Hp) 
p{p) 

(r, s) 
r(f) 

s h (f) 



A + P(0)/P(m) 

/ The a priori probability \ 
I that the signal is absent / 

The a priori probability" 
that the signal m(ii) 
(for undetermined //) 
is present 

d N nHv), 

P(p) exp [2 (g, m(ii)) - (m(fi),m(fi))] 
/The a priori probability \ 

density that the signal , 
^parameterization is 11. J 

3 ; 5,(1/1) ' 

dte 2nift r(t), 

One-sided detector noise 
power spectral density (PSD) 



(2.2) 

(2.3a) 

(2.3b) 

(2.3c) 
(2.3d) 

(2.3e) 

(2.3f) 
(2-3g) 
(2.3h) 



The likelihood ratio A (eqn. |2.3c|) is proportional to the a posteriori probability that 
a signal is present in the observed g(t). When that probability exceeds a given threshold 
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we assume a signal is present and try to determine jl. We characterize an observation of 
g(t) in terms of the mode (i.e., the maximum) of the distribution P(fi\g), denoted ft. The 
mode of P(/J,\g) is also the mode of the odds ratio A(/x). In terms of /2 the observation's 
signal-to-noise ratio (SNR) p 2 is! 



p 2 = 2{m(p),m(p)) . 



(2.4) 



Assuming that A(^t) is smooth and that its global maximum is a local extremum, ft satisfies 



m{/j,) — mln), 



dm 
dpj 



d In p , _ . 



- 2 /n,l^) 



(2.5) 



We assume that the noise n(t) is a normal random variable with zero mean; consequently, 



so are each of the (n, dm / dpj ) on the righthand side of equation [2.5| . Denote these random 
variables Uf. 



V; 




(2.6) 



Since the V{ are normal, their joint distribution is a multivariate Gaussian, characterized by 
the means VI, which vanish, and the quadratic moments 



ViVj 




\ a pi dpj 



(2.7a) 
(2.7b) 



Here we have used an overbar to indicate an average over all instances of the noise n(t). 
Since we have assumed that the noise is normal these averages can be evaluated using the 
ergodic theorem (c/. Finn [Q] for more details). In terms of Cjj (i.e., the inverse of C^ 1 ), 
the joint distribution of the Uj is given by the probability density 



p(u) 



cxp 



9 CijlSiV. 



(2tt 



M i,,il/2 



(2.8) 



This is also the joint distribution of the quantities that appear on the lefthand side of 
equation |2.5| . In Finn [B3] it was stated that for an observation characterized by ft, the 
probability volumes V(P) are given implicitly by 



2 Note that p 2 is quadratic in the signal strength. In the literature SNR is often used to refer to 
both p and p 2 . We avoid this ambiguity by using either p or p 2 as appropriate. 
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l-.l 



I /-\ dm \ dh\p 

2 I m{n) - m(», gj^W ) + "^tCA*) 

/ /-\ dm \ dliap 
2 ( m(» - m(/x), — (/Li) ) + — (fj,) 



(2.9) 



where 



P 



exp 



_ 9 CijiSiis. 



(27r) Jv det||Q 1 



1/2 



(2.10) 



and p, G V(P) with probability P (recall that P is conditional on the assumption that a 
signal is in fact present in the detector output). This is correct only when A (/a) has a single 
extremum, or near ft. In the most general case only a Monte Carlo analysis can determine 



the probability volumes V(P). No other results of Finn are affected by this correction 



C. The strong signal approximation 



The volumes V(P) are representations of a cumulative probability distribution function. 
Denote the corresponding probability density by P(Sfj,\fL): 



P(6»\fi) = 



/The conditional probability \ 
density that m is characterized 
by /2 + 5/x, given that A(/z) 

\has mode fi ) 



(2.11) 



In the limit of large p 2 , P(5[i\pL) becomes sharply peaked about p, and the determination 
of V(P) is greatly simplified. Suppose that p 2 is so large that for /2 e V(P) for all P of 
interest the difference m(Jj) — m(p) can be linearized in 5/j,, where 



5fi = \x — \x. 



We then obtain in place of equation 1275 



dm \ dlnp 
n, —(a*) ) - -7T— (A*) 



' dp 



.1 



dfji 



E^C^ = -2 

The random variables are related to the u by a linear transformation, 

consequently, the 5/j, are normal with means 

5 fit 



(2.12) 



(2.13) 



dlnp 



(2.14) 



^ Q \\\ p 



h3 



dfij 



(2.15) 
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and quadratic moments 



(5 ^ - 5 /J*) {6 m - S/ij) = dj. (2-16) 
The probability distribution P{5^\p) is thus a multivariate Gaussian (cf. eqn. |2.8| ): 



exp 



\ Eij <V (fyi - <^i) - fy; 



(27r) Ar det||C 



1/2 



(2.17) 



Note that the matrix Cy has now acquired a physical meaning: in particular, we see that 
the variances of of the Sfii are 



Of = - 8/1, 

= C 

and the correlation coefficients are given by 



f2.18l 



(TiUj 



[2.19) 



In this sense we say that C%j is the covariance matrix of the random variables Sfi. 

In the strong signal approximation, the surfaces bounding the volume V(P) are ellipsoids 
defined by the equation 



(2.20) 



'•j 



where the constant K 2 is related to P by 



d N x 



V. c' ■.<-.,■■■ i<- 

t—ii,] 1-3 — 



exp 



iy. . f^rWi 

2 <-—<i,j %3 



(27r) 7V det||Ci 



1/2 



[2.21] 



Finally we come to the question of when the linearization in equation |2.13| is a reasonable 
approximation. Three considerations enter here: 

1. The probability contours of interest {e.g., 90%) must not involve Sfi so large that the 
linearization of m(fi) — m(fi) is a poor approximation; 

2. The probability contours of interest must not involve Sfi so large that for Sfi G V(P), 
A(Sfi + fi) has more than one extremum or inflection point; and 

3. The condition number of the matrix C^ 1 must be sufficiently small that the inverse Cij 
is insensitive to the linearization approximation in the neighborhood of /i.i 



Recall that the relative error in 5a is the condition number times the relative error in Cr/: for 
a large condition number, small errors in Q- introduced by the linearization approximation can 
result in large errors in 5 a (cf. Golub & Van Loan p9[). 
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If the validity of the linearization procedure for a particular problem is doubtful owing to 
the violation of any of these conditions, then we must fall-back on either equations 2.5 , 2.9 
and [2.10| (if A(/x) has a single extremum), or a Monte Carlo analysis. 



III. PRECISION OF MEASUREMENT: 
BINARY INSPIRAL 

In this section we apply the techniques described in §|II| to the problem of measuring the 
characteristics of an inspiralling binary system in an interferometric gravitational radiation 
detector. We determine the gravitational radiation from an evolving binary system in the 
quadrupole approximation and evaluate the corresponding detector response. The amplitude 
of the response is a function of the (unknown and unknowable) relative orientation of the 
source and the detector, and we evaluate its mean square amplitude and a priori cumulative 
probability distribution. We also discuss the validity of describing binary evolution using 
the quadrupole approximation, arguing that while the result is certainly not good enough 
for use in actual data analysis, it is sufficient for the purpose of exploring the precision 
with which we will ultimately be able to characterize a binary system through gravitational 
radiation observations in LIGO/VIRGO-like interferometers. 

The general interferometer response to the gravitational radiation from an inspiralling 
binary system is a sinusoid with slowly varying amplitude and frequency. Using the station- 
ary phase approximation, we obtain an analytic expression for the Fourier transform of the 
response. We find that the SNR p 2 and the covariance matrix Cij can be expressed simply 
in terms of several moments of 5'^ 1 (/), the inverse of the interferometer noise PSD. 



A. Parameterization of the radiation waveform 

In order to express the response of a single interferometer to the gravitational radiation 
from a binary system (or any source), we define two coordinate systems: the source co- 
ordinate system and the radiation coordinate system. The binary system is most simply 
described in the source coordinate system. The ef-axis of this coordinate system is along the 
binary system's angular momentum. We choose the axes ef and of the source coordinate 
system to make the expression of the radiation directed toward the interferometer simple: 
the axis is chosen so that the unit vector n/ in the direction of the interferometer is in 
the A ef plane and in the positive ef direction (if the interferometer is along the polar 
axis, then there is no preferred direction for the ef axis). 

The radiation coordinate system has its ef axis in the direction n/, and its and 
axes are projections of the ef and coordinate axes normal to n/: 

e* = e s y , (3.1a) 
ef = x m. (3.1b) 

In the radiation coordinate system, the radiation propagating toward the interferometer is 
described by 
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h = h + e + + h x e > 



(3.2) 



where 



ef ® e.f + e J 



(3.3a) 
(3.3b) 



The quadrupole waveform 



In the quadrupole approximation^ 

M 



where 



h. 



2'— (1 + cos 2 ij (tt7W/) 2/3 cos ($ + *) 



h x = 4^ cos i (nMf) 2/3 sin ($ + *) 



cosi 



^Luminosity distance to binary^ , 
nj • ef , 

/Inclination angle of orbital 
i = angular momentum to line of 
\sight toward the interferometer. 

M. = (Chirp mass) , 
= (1 + z) fi 3 / 5 M 2 / 5 , 
T-t\ 5 ' 8 
5M ) ' 

i a$ 

2tt ' 



./ 



5 M 



1 3/8 



256T-t 



(3.4a) 
(3.4b) 



(3.5a) 
(3.5b) 

(3.5c) 



(3.5d) 
(3.5e) 

(3.5f) 



^ is the phase of the binary system at t = T, T is the Newtonian "moment of coalescence," 
z is the cosmological redshift of the binary system, and M and fi are the binary system's 
total and reduced mass P,[40f . The gravitational radiation frequency / is twice the systems 
orbital frequency. 



1 Here and henceforth we adopt units where G = c = 1 
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2. The interferometer response 



The response of an interferometer to the gravitational radiation field is a linear combi- 
nation of h. and h*: 



m 



F + h:+F x h y 



(3.6) 



where the antenna patterns F + and F x depend on the orientation of the interferometer with 
respect to the binary system EDI. We make no assumptions regarding the relative orientation 
of the interferometer and the binary; consequently, the general interferometer response is 

m(t; A,M,ip,T) = AM (iifM) 2/3 cos ($ + if)) , (3.7) 



where if) is a constant (distinct from \1/ in equations [3.4a| and |3.4b|) , and $ and / are given 
above in equations |3.5ej and |3.5f| . We will return in § [111 E| to discuss how A depends on the 
orientation angles through F + and F x . 



3. Radiation reaction and the quadrupole approximation 

LIGO/VIRGO-like interferometers are most sensitive to gravitational radiation with fre- 
quencies in the range 30-1000 Hz. This corresponds to binary orbital frequencies of 15- 
500 Hz. In this regime gravitational radiation reaction is the most important factor in 
determining the evolution of a binary system's orbit [31,[32]. In § 111 A 1 , where we describe 



the evolution of the binary systems orbit, only backreaction owing to the leading-order 
quadrupole radiation is taken into account This approximation neglects higher order 

effects (in both v/c and M/r) that contribute to the gravitational radiation luminosity and, 
consequently, the evolution of /. This has serious ramifications for the construction of the 



model detector response [41 



The detector response to binary inspiral is a sinusoid of slowly varying amplitude and 
frequency. The determination of the characteristics of the binary system is equivalent to 
finding the "template" response that is most closely correlated to the detector response. If 
the phase of the template drifts from that of the signal by as little as 7r radians over the course 
of the observations, then the correlation will be insignificant. Neutron star binary inspiral 
observations in LIGO/VIRGO-like interferometers will last for on order 2tt x 10 4 radians 
in phase; consequently, the phase advance can be determined to better than 1 part in 10 4 . 
The errors we have made in our template m(fx) (eqn. |3.7| ) by neglecting the post-Newtonian 
contributions to the evolution of the binary system lead to phase differences significantly 
greater than 2ir radians over 10 4 cycles observed. Consequently, the waveform model used in 
actual data analysis must be more accurate than that given by the quadrupole approximation 



5 Cutler, Finn, Poisson, and Sussman 42] have shown that successive post-Newtonian approxi- 
mations to the evolution and waveform converge very slowly upon the fully relativistic solution. 
Consequently, it may require an impractically high order post-Newtonian expansion to predict 
correctly the advance of the phase over the course of the LIGO/VIRGO observations. 
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Nevertheless, for the particular purpose of exploring our ability to detect and characterize 
binary systems by their gravitational radiation signature, we believe the quadrupole approx- 
imation waveform is a useful substitute for a more accurate waveform. The determination of 
the anticipated sensitivity of a LIGO/VIRGO-like interferometer to binary inspiral depends 
on the p 2 and Cij, and we expect that the quadrupole approximation provides good lower 
bounds on these. The predominant observable effect of the inclusion of post-Newtonian cor- 
rections is to change the rate that the wave frequency / advances, thus changing the elapsed 
phase of the wave over the period of the observation. The SNR p 2 depends approximately 
on the elapsed phase and the corrections, while large compared to 2tt, change this by only 
a small fraction of the total. Similarly, in the quadrupole approximation the rate at which 
/ advances depends exclusively on M. (c/. eqn. |3.5f ); consequently, to the extent that the 
corrections to the quadrupole formula depend on characteristics of the binary other than 
Ai (e.g., component masses and spins) these characteristics are observable and affect the 
precision with which M. can be measured. In this way we expect that the use of a more 
accurate waveform in our analysis will increase the estimated cr_A4 but have little effect on 
our estimates of 04, &ip, and <jt made using the quadrupole waveform and corresponding 
binary evolution. 



B. The stationary phase approximation 



In order to evaluate the SNR p 2 (eqn. |2.4j ) and the covariance matrix Cy (eqn. |2.7b| ) we 
must find the Fourier transform m of m (eqn. |3.7| ). We approximate m using the method of 
stationary phase. Given a real function of the form 



k(t) = A(t) cos$(t), 



(3.8) 



where <9$ jdt is a monotonically increasing function of t, the stationary phase approximation 
to the Fourier transform k(f) is 



*(/) 



where 



dtk(t)e 2nift 



\a [r{f)\ 

x exp 
[ **(-/) 



1/2 



2vr/!![T(/)] 

\ (27r/T(/)-$(T(/))+7r/4 



for / > 
for / < 0, 



(3.9a) 



(3.9b) 



T 



1 0$ . . 

T,m {t> 



t. 



(3.10) 



The validity of the approximation rests on the assumption that the amplitude A and the 
angular frequency d$/dt change slowly over a period: 
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d In A/dt 

d$/dt 
d 2 $/dt 2 

(d$/dtf 



< 1 and 

< 1. 



(3.11a) 
(3.11b) 



For the interferometer response m given in equation |3Ty] these two relations are equivalent 
to the single condition 



T-t /3tt\ 8/5 /256\ 3/5 



M 



> 



) (?)' 



14, 



(3.12) 



or, alternatively, 



5 8 \ 3/5 1 
~ 37' 



-(—-) 

7i V256 3tt/ 



For binary systems that will be observable by LIGO/VIRGO-like interferometers, 



fM = 6 x 10" 4 ^ 



M 1 

< 



100 Hz 1.2 M 37' 



(3.13) 



(3.14) 



consequently, the stationary phase approximation is a good one for our purposes. We thus 
have 



m(f) 



AM 2 

x 



384 



1/2 



nfMy 7/6 

exp{z [2nfT + ^ (tt/7W)~ 5/3 - ^ + f ] } for / > 0, ( 3 - 15 ) 

m*(—f) for / < 0. 



C. The signal-to-noise ratio 

Now suppose that we have analyzed the output g(t) of an interferometer for the signal 
m and found that the likelihood function is maximized for the parameterization {^4, Ai, ip, 
T}. The SNR p 2 is then given by (cf. eqn. 



p 2 = 2 (m{p),m{p)) 

= d^ A ^' 3f ^ (3 - 16) 

where, if we assume that we have access to g(t) for all t, 

/t/3 = / 4f / 7/3 ^(/) • (3-17) 



In practice, of course, data analysis is limited to a finite length sample of the interferome- 
ter output. During this limited interval, the signal "frequency" / (eqn. |3.5i|) of an inspiralling 
binary system ranges from f) to fh- We assume that the minimum of Sh(f) occurs in this 
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interval and not too near the endpoints. Then, just as the slowly varying amplitude and 
frequency of m permitted us to estimate m(f) using the stationary phase approximation, so 
it allows us to approximate p 2 for a finite duration observation by replacing the lower and 
upper limits in the integral expression for / 7 / 3 (eqn. |3.18|) by f and fh- The lower bound 
fi is determined by the duration of the data stream being analyzed while the upper bound 
fh is determined by the coalescence of the binary components. We further assume that the 
ft/3 does not change significantly as fh — > oo; thus, we have 



/7/3 



df \f 7/3 S h (f) 



(3.18) 



We discuss our choice of low-frequency cut-off f in §[IV A 6| and again in §VD 



D. The covariance matrix 



Turn now to the calculation of the covariance matrix Cij (cf. eqn. |2.7b| ). Instead of a 
parameterization in terms of A and Ai, it is more convenient to introduce rj and ( defined 

by 



A(l + rj) = A 
M(1 + C) = M, 



(3.19a) 
(3.19b) 



where A and M. are the modes of the observed distribution of A and A4. 

Given our expression for m(f), we can evaluate all of the elements of the symmetric 
matrix in terms of the frequency moments fa defined by 



Sa = f- 



7/3 



df \fSh(f) 



-i 



(3.20) 



To express dj, only f* for (3 G {17/3, 4, 3, 7/3, 1/3} are needed. In terms of these moments, 



/ 1 



C" 1 



25 
36 



1 + 



/l7/3 



4096 



10/3 



V 



J.2S 



1 



573 





5tt fz 



\ 



64 / -~-\S73 
( ttM) 

-27t/4 /3 

(2vr) 2 /i/ 3 J 



(3.21) 



where the indices are ordered r], (, ip, and T. We have inverted equation |3.21| to find 
and so the variances and correlation coefficients describing the distributions of 77, (, ip, and 
T (cf. eqns. |2.18| and |2.19p. The variances are 



a 



A 



A 2 



4096 / n \io/3 /- - 2 x 



-2 



(3.22a) 
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M 2 
16384 



25Ap 2 

2 _ /l7/3/l/3 ~ ^3 

2 = 7l7/3 ~ 7 4 
47T 2 Ap 2 

and the correlation coefficients are 



(/1/3 - fl/3) 



10/3 



64 (f\, 3 - / 1/3 ) (ttM) 



10/3 



9 A + 4096 (ttM) 10/3 (/ 1/3 - /J /3 ) (/ 1/3 - /J /s ) (ttA<) 

5/3 



10/31 V2' 



64 ( f 1/3 f 4 - f 4/3 f 3 (ttM 



{ [9 A + 4096 (vr^) 1073 (/ 1/3 - /J /3 )] (/ 17/3 / 1/3 - / 3 )} 

6 1(7,7;, 7,) (■-m' : 



1/2 • 



9 A + 4096 (vr^) 1073 (/ 1/3 - /* /3 )] (/ 17/3 - /5} 1/2 



r Ci> — r Mtp 



(74/3/3 /1/3/4) 



where 



(/1/3 ^4/3) (/17/3/1/3 / 

_ /3 ~~ f if 4/3 

{fl/3 ~ ^4/3) (/l7/3 ~~ ^4) 
/17/3/4/3 ~~ ^4/3 



1/2 ' 



(/l7/3/l/3 ^3) (/l7/3 A) 

a = (7i /3 - 7! /3 ) 7it/3 - 7i /3 74 + (27 4 74/ 3 - 7s) 7 3 > o. 



(3.22b) 
(3.22c) 
(3.22d) 



(3.23a) 



(3.23b) 



(3.23c) 



(3.23d) 

(3.23e) 
(3.23f) 

(3.24) 



E. Properties of „4 

Despite the fact that we cannot measure the orientation angles relating the interferometer 
to the source, we still need them in order to assess the interferometer's sensitivity. Recalling 
that 
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m = F + h + + F x h x , 
we express F + and F x according to the following convention: 



(3.25) 



1. Assume that the interferometer arms are the same length and that they meet in right 
angles. 

2. Define a right-handed coordinate system with one interferometer arm along the x-axis 
and the other along the y-axis. Denote the unit vector in the direction of the x arm 
by 1 and the unit vector in the direction of the y arm by m. 

3. Let the position of a source in the sky be given by the polar angle 9 and the azimuthal 
angle 0, and denote the unit vector pointing toward the source by ris {i.e., — nj). 

4. The interferometer responds linearly to the radiation field h, so its response m can be 
represented by a tensor R such that 



m = R:h. 



For our interferometric detector 



R=-(l<g>l-m<g>m). 



(3.26) 



(3.27) 



5. Assume that axis makes an angle ( with the axis 1, i.e., 

= 1 cos C + m sin (. 



(3.28) 



With these conventions, the antenna patterns are given by |g0 

F + = R:e+ 

- cos 2( (l + cos 2 cos 2<p — sin 2£ cos 6 sin 



F x = R:e x 



sin 2((l + cos 2 9) cos 20 + cos 2£ cos 9 sin ! 



and the amplitude A 2 may be written 



A 2 



dr 2 



F'i (1 + cos^ iY + 4F 2 cos 2 i 



It is convenient to denote the angular dependence of A by 6 : 

, 2 



e 2 = 4 



F'f (1 + cos^r +4F 2 cos 2 i 



(3.29a) 
(3.29b) 

(3.30) 

(3.31) 



The range of 6 2 is < 6 2 < 16. 
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The SNR p 2 and covariance matrix dj (eqns. |2]4] and [2.7b| ) both depend on A 2 , which 



is in turn a function of the (unknown) relative orientation of the source and interferometer 
through O 2 . In order to evaluate the expected p 2 (or the expected Cij) of a source at a given 
distance we need to know some properties of the probability distribution of 2 . 

Since 6 2 depends on the angles 9, 4>, i, and ( (cf. eqn. |3.31| ), the a priori distribution of 
2 depends on the a priori distribution of these angles. These distributions are all known: 
in particular, cos# and cosz are uniformly distributed over the range [—1, 1] and <ft an d C 
are uniformly distributed over the range [0, 2tt). Making use of the definitions of F + and F x 
(eqns. |3.29a| and |3.295| ), we find that the mean square of is 



64 

02 = 25" (3 ' 32) 

The distribution of 2 is not symmetric, however: in fact, its mode is zero and larger values 
of 2 are much less likely to occur than smaller ones. We have determined the cumulative 
distribution function of 2 using a Monte Carlo analysis; we give the percentiles of the 
distribution in table |. In performing these calculations, we used Knuth's portable random 
number generator [J3| as implemented by Press et al. JPJ (i.e., their RAN3). The results 



in table [I] are based on a sample of 10 7 points in the {cos^, cosi, 4>, (} parameter space. 
The corresponding values of 2 were sorted into bins and the reported percentiles are the 
rounded bin centers. 

Note from table [j] that significantly more than half (i.e., approximately 65%) of the 
inspir ailing binary systems will have 2 less than 2 . The skew of the distribution toward 
smaller 2 plays a significant role when we estimate the range of the interferometer (cf. 

gym 



IV. APPLICATION TO LIGO 

In § piT| we found expressions for the SNR, the variance, and correlation coefficients cor- 
responding to the detection of a signal characterized by A, A4, x/j, and T . In this section we 
join those expressions with the design characteristics of LIGO / VIRGO-like interferometers 
to obtain estimates for the sensitivity of a realistic interferometric detector to inspiralling 
binaries. 



A. Noise and the LIGO interferometers 

The characteristic of the interferometers enter our analysis solely through the strain noise 
PSD Sh(f)- The dominant contributions to Sh(f) for the interferometer configurations that 
will be used to search for inspiralling binaries are from seismic, thermal, and photon shot 
noise. In this subsection we summarize these contributions to the overall noise PSD that we 
use in our calculations. 
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1. Photon shot noise 



At high frequencies Sh(f) is dominated by photon shot noise. The shot noise depends 
on the interferometer arm lengths, laser power and wavelength, mirror reflectivities and 
configuration in a complicated fashion. In all of our calculations we have assumed that 
LIGO/VIRGO will be equipped with Fabrey-Perot cavity interferometers, and we have used 
the analysis of Krolak, Lobo, and Meers |28| , cf. their eqns. 2.11-22] to describe the photon 
shot noise. This analysis is general enough to encompass non-recycling, standard recycling, 
and dual recycling interferometers. Except in § |1V B| we always assume that the instrumen- 
tation makes use of standard recycling techniques. For standard recycling, the photon shot 
noise is given approximately by Thorne j|0], eqn. 117c] or Krolak, Lobo, and Meers f28| , 
eqn. 3.5] 



oshot 



(/) = ^ 




(4.1) 



where A 2 describes the mirror losses, Jo is the laser power, r\ is the quantum efficiency of 
the photo-detector, A is the laser wavelength, L is the length of interferometer arms, and f c 
is the recycling "knee" frequency. 

The simplest way in which the observer can change the noise characteristics of the 
LIGO/VIRGO instrumentation is by changing the recycling knee frequency f c . The gen- 
eral trend is that, as the recycling frequency increases the bandwidth increases while the 
sensitivity across the bandwidth decreases [40], fig. 9.13]. 



2. Thermal noise 

At lower frequencies, off-resonance thermal excitations of the test mass suspensions and 
internal modes of the pendulum masses either dominate or provide important contributions 
to the noise. We approximate the suspension noise by focusing only on the pendulum 
mode (ignoring both torsional and violin modes). If the dissipative force in the pendulum 
suspension is due to friction, then the strain PSD of a single test mass m with resonant 
frequency fo and quality factor Q at temperature T is given by 



gpend, ,x = k B Tf 

' " ' ' 2~>m(.) (l I 2 \(P - f 2 ) 2 + (/ fo/Qo) 21 



(note that we use T for both the temperature and the Newtonian "moment of coalescence;" 
nevertheless, the meaning of the symbol in any given context should be clear). Each arm of 
the interferometer has a pendulum degree of freedom at each end and the noise from each 
degree of freedom is independent; thus, the total noise PSD is a factor of four greater than 
this (note that Dhurandhar, Krolak and Lobo [22, eqn. 2.10] neglect this factor of four in 
their analysis and have several typographical errors in their formulae). The total thermal 
suspension noise PSD is given by 

sr p (/) = 4sr nd (/)- (4.3) 
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The primary dissipative force acting on the pendulum may not be friction, however: it 
has been suggested [^J that the dissipation is due instead to a phase lag between the stress 
and the strain in the pendulum suspension. If this is the case, then the noise PSD is different 
than that given above. This is not an issue for the initial interferometers, but will be for the 
advanced ones [HJ . The nature of the dissipation is far from settled, and in the absence of 
a consensus we have used the form given in equation fi~3 . 

Off-resonance thermal excitations of the vibrational modes of the test masses will also 
be a significant source of noise in the LIGO/ VIRGO interferometers. Here we consider only 
the fundamental vibrational mode of each test mass. The contribution to the noise PSD 
has the same form as the thermal suspension noise S*™ sp (/) (and is subject to the same 
controversy), only now the resonant frequency and oscillator quality that enter are those of 
the test masses: 

2k B Tf [nt 



ST(f) 



n*mQ int L* (P - fl t f + (ff int /Q, 



(4.4) 



3. Seismic noise 



Seismic noise will dominate Sh(f) at low frequencies. Saulson [i5| has surveyed the liter- 
ature on the seismic displacement noise PSD S x (f) and finds that it is roughly proportional 
to /~ 4 in the range 1/10 Hz < / < 10 Hz. Consequently, if the LIGO/VIRGO test mass 
suspensions were coupled directly to the Earth, then 

Sf smic (f) = , (4.5) 

(/<? - Pf + (fh/Qo) 2 

where fo is the pendulum mode frequency, Qo is pendulum quality, and S' Q is a proportion- 
ality constant with units of Hz 7 . Note that at frequencies above the pendulum frequency f 
the seismic strain noise is proportional to /~ 8 while below fo is is proportional to f~ 4 . 

In the actual LIGO/VIRGO interferometers, the pendulum suspensions will be isolated 
from the Earth by a mechanical circuit that is a series of several highly damped oscillators 
4"6| , |4"?|| . Each oscillator in this series circuit will introduce four poles in the strain noise PSD 



near the pendulum frequency; consequently, the actual seismic noise contribution will be 
much steeper than /~ 8 at frequencies greater than f and much more complicated near the 
resonant frequencies of the mechanical circuit. It is proposed that seismic isolation in the 
initial interferometers be provided by a five stage circuit (where the final stage of isolation is 
the pendulum suspension itself) [[|^, §V 3 b and appendix D]; consequently, for frequencies 
much greater than f the seismic strain noise PSD is expected to be proportional to f~ 24 . 
Nearer the resonant frequencies of the isolation circuit the dependence is quite complicated 
as the poles in the response function are not at zero frequency, but at complex frequencies 
with real parts near 1 Hz. The strain noise PSD below the resonant frequencies of the 
isolation circuit remains proportional to the displacement noise PSD S x (f). 

In our calculations we use the following crude estimate for the seismic strain noise PSD: 

sr mi v) = (/ f J ^ )10 - (4.6) 
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The proportionality constant So has units of Hz 23 . This estimate scales correctly with 
frequency above and below /o, though it fails near / . This failure is unimportant since 
(except in §V D| ) we assume that fi > f (the choice of cut-off fi is discussed in §1V A 6| ). 

The amplitude of the noise (as reflected by the proportionality constant S ) depends 
on the detailed nature of the seismic isolation circuit and the properties of the seismic 
displacement noise at the interferometer site. The LIGO design goals for the initial and 
advanced interferometers are 

^seismic ( 4Q = 5 susp( 4Q R ^ + ^int^Q R ^ interferometers ) ( 4 . 7a ) 

5 seismic( 10 Hz ) = S ^p^q Hz ) + 5 mt^ (advanced interferometers) . (4.7b) 

We fix So by these relationships. Where we discuss the advanced LIGO detectors, we use 
the same condition on the seismic noise as for the advanced interferometers. 



4- Quantum noise 

In addition to the primary noise sources discussed above, we have also included a con- 
tribution whose origin is quantum mechanical and rooted in the Heisenberg uncertainty 
principle. When we observe a signal of frequency / in an interferometer, we are measuring 
the periodic motion of the end-masses at that frequency. Since the motion is periodic, this 
is equivalent to a simultaneous measurement of the momentum and localization of the end 
masses, and the precision with which we can make this measurement is subject to the usual 
quantum mechanical limits. In our calculations, we use the form of the quantum noise given 



by Thorne |40| , eqn. 121]: 



sr nt (/) = *% T2 , (4.8) 



where m is the mass of the LIGO pendulum bobs. 



5. Noise source summary 

In table |I| we give the instrument characteristics we have assumed in our calculations. 
Two sets of values are given, corresponding to estimates for LIGO initial and advanced 
instrumentation. These estimates have been culled from the literature |]|§, the LIGO 
proposal fj6 |, and personal communication with members of the LIGO project [|l^,Q. In 



terms of the noise sources discussed in the previous subsections, the noise PSD Sh(f) we 
use in our calculations is 

s h w = st t (f) + s^(f) + sr p (f) 

+ Sf smic (f) + Sr ant (f). (4.9) 

As a companion to table [Tl] and as a graphical illustration of how all of the noise sources 
discussed above act in concert to determine an interferometer's noise characteristics, we show 
our approximation to the anticipated Sh{f) for both the initial (fig. |1|) and advanced (fig. 0) 
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instrumentation. The contributions to Sh(f) from each of the influences discussed above are 
shown as dashed lines and Sh{f) is shown as a solid line. Both figures show interferometers 
configured to operate in standard recycling mode. In figure [1], corresponding to the initial 
interferometers, the recycling frequency f c is 300 Hz, while in figure |2| (corresponding to the 
advanced interferometer design) it is 100 Hz. 



6. Choosing the low frequency cut-off fi 

The elements of the covariance matrix dj depend on the moments f 7 / 3 , f 17 / 3 , / 4 , / 3 , / 4 / 3 , 
and fy 3 . In turn, these depend on Sh(f) and the low frequency cut-off fi (see eqns. |3.18 



and p.20|) . Our calculations of these moments have assumed a low frequency cut-off of 10 Hz, 
corresponding to the last several minutes in the inspiral of a binary neutron star system. In 
table m we give the frequency moments f 7 / 3 , fn/ 3 , / 4 , f 3 , / 4 / 3 , and f 1 / 3 for two cases of 
interest . 



7. The LIGO detector 



The two LIGO interferometers, though separated by several thousand miles, share nearly 
the same orientation in space: the planes defined by the detector arms are nearly parallel, 
and the arms themselves are nearly parallel. Consequently the network acts like a single 
interferometer of greater sensitivity than either of its components. If the noise in the two 
component interferometers of the LIGO detector is uncorrelated and described by Sj£ (/), 
then the effective PSD of the more sensitive single interferometer is Sh(f) = S^\f)/2. 
Consequently, the effective PSD for the LIGO detector in the limit that the interferometers 
share the same orientation is also given by figure |2], but with the scale reduced by a factor 
of 2 -1 / 2 . In making this approximation we are ignoring the differences in arrival time of 
the gravitational radiation signal at the two interferometers. In the following sections when 
we refer to the LIGO detector (as opposed to a LIGO/VIRGO-like interferometer) we are 
actually referring to a single interferometer whose noise PSD is 1/2 that of the advanced 
interferometer design. 



B. Signal-to- noise ratio 



As discussed in §11 B, we decide whether or not a signal is present in the output of the 



detector by comparing the likelihood ratio A to a pre-determined threshold. In this regard, 
the SNR p 2 is an acceptable surrogate for A; i.e., we can choose a threshold p^ (which 
may be a function of fi) to compare with p 2 . Then, if p 2 > p\ we assert the presence of a 
signal while if p 2 < p\ we conclude that the detector output is only noise. The choice of 
threshold is a delicate matter: on the one hand we want a high threshold to minimize the 
probability that we misidentify noise as signal; on the other hand, we want a low threshold 
to minimize the probability that we misidentify a real signal as noise. We will consider the 
proper choice of the threshold p\ in a later paper; now, however, we assume only that the 
threshold depends weakly on the detection strategy and /x. 
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The amplitude SNR p may be expressed 



M 



5/6 



i.2M ; \e 50% ) 

( (17.0Mpc/d L ) 
x I (308. Mpc/d L ) 
[ (436. Mpc/d L ) 



'7/3 



initial interferometers 
advanced interferometers 
advanced LIGO detector 



(4.10) 



where 



F7/3 (fc, fi) 



A/3(/o fi, 



.1/2 



/ 7/3 (100Hz,10Hz)_ 
' (/7/3 / 5.331 x 10 42 Hz- 1/3Xl 1 

(/ 7/3 /1.747 x 10 45 Hz 
^ (A /3 /3.493 x 10 45 Hz 



-1/3 
-1/3 



1/2 
1/2 



initial interferometers 
advanced interferometers 
advanced LIGO detector 



(4.11a) 



(4.11b) 



In figure [| we show F 7 / 3 for both the initial and advanced LIGO instrumentation. Two curves 
are shown: one for the initial interferometers and one for the advanced interferometers. Each 
curve assumes a low-frequency cutoff fi of 10 Hz (corresponding to the last several minutes 
of binary neutron star inspiral). For the advanced instrumentation the detection strategy 
that maximizes p has f c = 100 Hz (where F 7 / 3 = 1), while for the initial instrumentation 
f c = 300 Hz (where F 7 / 3 = 1.3). These correspond to the choice of f c in figures [I] and |2| 
showing the detailed breakdown of Sh(f) for the initial and advanced interferometers. 

Now consider the case of resonant dual recycling |T9|j2T|j24| . The photon shot noise in 
an interferometer operating in a resonant dual recycling mode is proportional to (cf. |28 
eqn. 3.7 but note that their approximate expression has several errors; see also |i~9lj2"4jl ) 



(/)« 



7 - fn 



(4.12) 



where Af and f n depend on the reflectivities of certain mirrors in the experimental appara- 
tus. This shot noise PSD is large except in a narrow band about f n where it is very small. 
The bandwidth Af and the central frequency /„ of this "notch" can be adjusted, and the 
size of S^ ual in and out of the notch will vary depending on the / and Af. In comparison, an 
interferometer operating in standard recycling mode has a nearly constant shot noise PSD 
for frequencies below the knee frequency f c , with the noise PSD increasing as f 2 for frequen- 
cies greater that f c (cf. eqn. [O] and figs. [I] and H). For \f — f n \ ~ Af, dual recycling cuts 
a notch in S^ ual . Elsewhere, Sf^ ot is lower than that for dual recycling (assuming / ~ f c ). 
On the basis of numerical investigations of the interferometer response and the waveform of 
inspiralling binary systems, Krolak, Lobo, and Meers |^SJ suggested that dual recycling with 
/„ = 100 Hz and Af ~ 6 Hz is superior to standard recycling for the observation of binary 
systems. Their analysis assumed that Sh{f) is infinite below 100 Hz and that only photon 
shot noise is important above 100 Hz: in particular, they did not consider the noise owing 
to the thermal excitations of the test mass vibrational modes (cf. § |LV A2| ). When S^ t (f) 
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is included in the analysis, then the advantage of dual recycling in this regime is lost: the 
notch is "filled in" by the thermal noise almost to the level of standard recycling photon 
shot noise (cf. their fig. 5 and our figs. |I]or 0), and everywhere else the noise is much greater 
than is the case for an interferometer operated in standard recycling mode. 

Dual recycling may still be useful at much higher frequencies where the photon shot noise 
is a standard recycling configuration is much larger than the thermal noise (e.g., 1.5 KHz). 
Observations at such high frequencies may prove useful for detecting the actual coalescence 
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V. ASTROPHYSICAL IMPLICATIONS 
A. Source rate 



Chernoff and Finn fl£| have shown that the observed differential rate dN /dM.dA of 



inspiralling binary systems depends on the cosmological model; consequently, it can be used 
to determine the Hubble constant H, the deceleration parameter q, and otherwise distinguish 
between cosmologies. Here we are interested in the total rate of observed binary inspiral as an 
estimate of the sensitivity of an interferometer, and in this subsection we estimate that rate 
ignoring cosmological effects. We refer the reader interested in a rate calculation consistent 
with an expanding universe and taking into account evolution of the binary population and 
distribution of j\A in binaries to Chernoff and Finn [[Tjj]. 

Assume that the rate density (number per unit co-moving cosmological volume per unit 
time) of inspiralling neutron-star binary systems is a constant Af and that the variation in 
neutron star masses is small so that A4 is approximately equal to 1.2 M Q (corresponding to 
two 1.4M Q neutron stars). The expected total rate N of systems whose SNR p 2 is greater 
than Pq is 



N = / drAitNr 2 P p 2 (r) > pi 
Jo L 

r°° • ( v 2 \ 
= / drA-nMr 2 P 6 2 > — 
Jo \ r n J 



AnjVrlJ dxx 2 P(Q 2 > x 2 ) (5.1a) 

fivn\ (5.ib) 



where 



P ^0 2 > x 2 ^j = (The probability that B 2 is greater than x 2 ) (5.2a) 

^7/3 



^0 = 



96vr 4 / 3 po J 
M \ 5/6 



1.2 M n 




24 



X 



13.0 Mpc initial interferometer 
236. Mpc advanced interferometer 
334. Mpc advanced LIGO detector, 



n = r 



dxx 2 P (e 



2 ^ 2 
> X 



1/3 



(5.2b) 
(5.2c) 



[recall that we have ignored cosmological effects (cf. Chernoff and Finn jllj) in eqn. |5.1a 
Using the cumulative distribution function for 2 (cf. §111 E| and tbl. |), we find that 



J™ dxx 2 P(& 2 > x 2 ) = 1.84; 



(5.3) 



hence, 



K 



M 
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7/3 



1.2M ; VPo, 

23.0 Mpc initial interferometer 

x < 417. Mpc advanced interferometer 

589. Mpc advanced LIGO detector. 



(5.4) 



Phinney [[| has given estimates for the number density of sources per unit time (Af) 
based on observational and theoretical arguments. These estimates range from an ultra- 
conservative 6 x 10~ 10 Mpc _3 yr _1 , to a conservative 8 x 10~ 8 Mpc _3 yr~ 1 , to an upper limit 



of 6xlO- 5 Mpc~ 3 yr _1 - They are based on the statistics of local populations of binary pulsars 
and type lb supernovae, and the large range reflects both the small size of the local sample, 
uncertainties in our understanding in the evolution of binary systems, and uncertainties in 
the selection effects at work in determining the fraction of the local systems we have direct 
knowledge of. If we take the typical threshold po to be 8, then we find that expected rate of 
detections of inspiralling binary systems is 

5/2 / n \ 3 



N 



Af 



M 



8 x 10- 8 Mpc _3 yr _1 
' 4.1 x 10~ 3 yr- 
x I 24. yr -1 
69. yr _1 



1.2M 




0/ 

initial interferometer 
advanced interferomter 
advanced LIGO detector. 



(5.5) 



As commented earlier, to maximize the rate at which binaries are detected we need to choose 
f c in order to maximize -F7/3. 

For a binary system consisting of a 10 M black hole and a neutron star, A4 ~ 3, and for 
a binary system consisting of two 10 M black holes, M. ~ 9 (recall that we are neglecting 
cosmological effects). Consequently, for these neutron-star/black-hole (black- hole/black- 
hole) binaries 7^go% — 2 Gpc (5 Gpc) for the advanced LIGO detector. The situation 
for determining the rate at which such systems will be detected is a bit more complicated. 
Phinney argues that black-hole/black-hole and black-hole/neutron-star binaries form at rates 
comparable to the neutron-star/neutron-star merger rate; however, the fraction which merge 
depends on the model dependent details that vary greatly M, so no reliable estimate of the 



coalescence rate is available for use with eqn. 5.5 
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B. Range 



An important measure of the sensitivity of a LIGO/VIRGO-like interferometer is its 
"range," i.e., the distance to which sources can be observed. The definition of the range is 
subtle. Not all inspiralling binaries within, e.g., a distance 1Z will be identified as such: for 
some, O will be less than 7Z/r , the corresponding SNR p 2 will be less than the threshold 
Pq, and the signal will be dismissed as noise. Similarly, not all inspiralling binaries outside 
a distance 7Z will fail to be identified by the interferometer: for some G will be greater than 
TZ/ro, the SNR p 2 will be greater than p 2 ,, and the signal will be identified as coming from 
a binary system. Since the range is a slippery concept, we define a range function 7Z^ such 
that a fraction 7 of the observable sources fall within the distance 7Z^: 

_ f^ /r °dxx 2 P(Q 2 >x 2 ) 
7 f£°dxx 2 P(Q 2 >x 2 ) ' 1 ' 



The quantity 

f*dxx 2 P{Q 2 >x 2 ) 



dxx 2 P (G 2 > x 2 ) 



(5.7) 



is tabulated in the third column of table |, and we show 7 as a function of IZ^/IZ in figure |j. 
Note that 

!37.2Mpc initial interferometer 
673. Mpc advanced interferometer 
952. Mpc advanced LIGO detector, 

(5.8) 

i.e., for the advanced LIGO detector approximately 7 sources per year will be observed whose 
distance is greater than 950, Mpc. Like the rate, the range is sensitive to the detection 
strategy so that if we wish to maximize the sensitivity of the interferometer to either we 
need to choose a detection strategy that maximizes F 7 / 3 . 



C. Standard deviation and correlation coefficients 

First consider the measurement of A. For all A4 and f c relevant for both the initial and 
advanced LIGO/VIRGO-like interferometers, the fractional standard deviation a v (cf. |3.22a|) 



of the waveform amplitude A is p x ; consequently, the maximum fractional one-sigma un- 
certainty in the determination of A is 

-4 = 0.125— 5.9 
A Po V ' 

for both the initial and advanced interferometers. Additionally, the correlation coefficients 
r_4j all have magnitude less than 10~ 4 , indicating that A is statistically independent of M., 
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ip and T (i.e., the random errors in measurements of A owing to detector noise are not 
correlated with the corresponding errors in the measurement of Ai, if), or T). 

Now turn to the measurement of M.. Before discussing the exact results obtained with 
equation |3.22b| , we give a heuristic derivation of the precision with which M. can be deter- 
mined. Recall that the phase $ of the gravitational wave signal is given by 



T-t 
, 5M 



5/8 



(5.10) 



The observation encompasses approximately the last ten minutes in the life of the binary 
system, during which time the phase advances by 



A$ = 7.4 x 1C (-^-H£^ rad. 
VlOm M 



The argument of the exponential in the odds ratio (eqn. (2.3d| ) is 
2 (n, rn(fi)) + 2 (m(£t), m(/x)) - (m(/x), m(/z)) . 



(5.11) 



(5.12) 



The contribution owing to the term 2 (m(/2), m(fi)) is much greater than that owing to the 
noise; consequently, to a good approximation the odds ratio will be maximized where this 
quantity is maximized. Ignoring the frequency dependence of Sh{f) the term (m(/i), m(/z)) is 
the correlation between two sinusoidal functions of the phase, and is large only as long as the 
advance in phase of m(/x) is within approximately it radians of the advance in phase of m(/2) 
over the course of the observation. Since A$ depends only on Ai, we have 7r > |A$ — A$| 
or 



5M < 4:71 




(5.13) 



where 5M is M — M . 

Return now to consider the exact results. Equation |3.22b| gives the fractional standard 
deviation aj^/Ai in terms of the frequency moments fa and p: 



M 



M 



5/3 



1.2 M 



2.08 x 10~ 
2.20 x 10" 



initial interferometer, 
advanced interferometer 



(5.14) 



where 



fl/3 fi/3\ I ( fl/3 f A/3 



A 



/ 



A 



/ c =100 Hz,/i=10Hz 



1/2 



(5.15) 
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Given a threshold po such that p > po for all observed sources, equations |5.14j and |5.15| give 
the maximum fractional standard deviation in the measurement Ai for any binary system 
observed with LIGO/VIRGO-like interferometers — a phenomenal precision. In interpreting 
eqn. |5.14| , note that a M / M. is inversely proportional to p, and recall that the SNR p of a 
binary system observed in an interferometer of the advanced design is approximately 26 
times greater than the SNR of the same binary observed in a detector of the initial design 
(cf. eqn. |4.10| ). The results for the advanced LIGO detectors are identical to those for 



the advanced interferometers, except that the amplitude SNR p for a binary observed in 
the advanced detector is 2 1//2 greater than that for the same binary observed in a single 
advanced interferometer. 

In figure [| we show for both the initial and advanced interferometers. For the ad- 
vanced interferometer the total variation of is approximately 20% as f c ranges from 50 Hz 
to 1 KHz, while for the initial interferometer the variation is approximately 15%. The op- 
timum recycling frequency for measurement of Ai is that which minimizes and we see 
that this is very different than the choice which maximizes the number of binaries observed 
(cf. § [i V B| and fig. |3]): in fact, the optimal choice of f c for the detection of binaries (in either 
the initial or advanced interferometers) is close to the worst possible choice of f c for the 
precise measurement of Ai. 

As we have pointed out, with observations in a single gravitational wave interferometer 
the location of the source on the sky cannot be determined. If, as has been suggested, some 
coalescing binaries result in 7-ray bursts |]49|j5[|, then burst observations may be used to 
localize the binary system in the sky. The identification between a gravitational wave burst 
from orbital decay (which takes place before actual coalescence of the binary components) 
and a 7-ray burst (which takes place at the time of coalescence) depends on the accuracy 
with which we can measure the time of arrival of the 7-ray burst and the "moment of 
coalescence" T: in all events, T will be within seconds of the actual moment of neutron star 
disruption and the emission of the 7-ray burst. Consequently, we need to know T, the rate 
of detected binary coalescence, and the rate of 7-ray bursts (the latter both assumed to be 
Poisson distributed in time) in order to evaluate the probability that a correlation in time 
between a 7-ray burst and a gravitational wave burst is coincidental. The accuracy with 
which we can determine T is given by o~t (cf- |3.22d| ): 



o~t 




1.54 x 10" 4 s 
3.00 x 10" 4 s 



initial interferometer, 
advanced interferometer, 



(5.16) 



where 



/17/3 [4 
A 



./ 



17/3 



A 



/ c =100 Hz, fx =10 Hz 



1/2 



(5.17) 



The results for the advanced LIGO detector are the same as those for a single advanced 
interferometer. The factor Ey varies by approximately a factor of 2.5 over the range 50 Hz < 
f c < 1000 Hz, and is shown (for both the initial and advanced interferometers) in figure M. 
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Again, the optimal interferometer configuration for precision measurements of T is very 
different than that for detection of inspiralling binaries. 

The parameter if> depends on the orientation of the source and the detector and the 
phase of the binary systems orbit at t = 0. For completeness, we also give the precision with 
which if) can be measured: 



X 



0.257 rads initial interferometer, 
0.338 rads advanced interferometer, 



(5.18) 



where 



/l7/3 ~ fi \ i ( f 



I 



17/3 



A 



/ c =100 Hz, fi =10 Hz 



1/2 



(5.19) 



The results for an advanced detector are the same as those for a single advanced interfer- 
ometer. 

The correlation coefficients rMip, r^T, and r^y are nearly independent of f c for both the 
initial and advanced LIGO-like interferometers. As mentioned above, the statistical error in 
A is essentially uncorrelated with that in A4, ip, or T (i.e., the correlation coefficients 
are for all < 10~ 4 ). Figure |7| shows the remaining correlations coefficients r M ^, r MT , and 
r^T for the initial and advanced interferometers. 



D. The low frequency cut-off 

In order to evaluate the covariance matrix we needed to compute the six frequency mo- 
ments /V/3, fi 7 / 3 , f 4 , f 3 , / 4 / 3 , and fy 3 . The evaluation of all of these is straightforward; 
however, the calculation of / 17 / 3 deserves special attention: at frequencies below the pendu- 
lum frequency of the LIGO masses the seismic noise PSD is proportional to f~ 4 (cf. eqn. |4.5|) : 
consequently, f 17 / 3 diverges as /; approaches (cf. eqn. |3.20| ). 

For any particular application we never encounter the divergence: there is always a low 
frequency cut-off in the integral |3.20| corresponding to the finite period of the observation. 
Even if we had access through an interferometer to the entire life history of a binary system, 
our model for its evolution is relevant during only a small part of its lifetime: for example, 
we have assumed that the orbit is circular for all times (when in fact gravitational radiation 
may be responsible for circularizing it), that the evolution of the orbit is due exclusively to 
the gravitational forces acting between the two components, and that the two bodies are 
bound in a binary into the infinite past. Similarly, our model of the detector noise PSD Sh(f) 
is not necessarily valid at very low frequencies. Nevertheless, the divergence of f\ 7 / 3 tells us 
something interesting about the observation of a binary system in gravitational radiation 
and it is worthwhile to spend a few moments understanding its origin. So, for the purpose 
of understanding this divergence we assume that we know Sh(f) as / tends to 0, that our 



29 



model for the waveform from a binary is correct throughout its life history, and that the 
lifetime of a binary extends into the infinite past. 

Turn first to the related example of a strictly monochromatic signal of frequency /o : 

m(t) = Acos2n f t (5.20) 

for all t. The noise PSD of the detector is given by Sh(f)- If we observe the signal for a 
finite period of time r, the SNR p 2 is 

No matter how small A 2 is compared to the power in the noise, the signal can always be 
discerned given a long enough observation time. The situation with binary inspiral is similar: 



consider the transformation of the time coordinate |25|] 



^-M(—) . (5.22) 

In terms of t', the signal can be expressed 

AM ( M\ 5/32 

h{t) = £L±(- J -±) cos($ + ^) (5.23a) 

$(t') = iTxfi! (5.23b) 

/' = — rr- (5.23c) 

Thus, the signal from an inspiralling binary is very much like that from a monochromatic 
source of radiation, save that 

1. The signal amplitude tends to zero as the stretched time t' tends to — oo; 

2. The detector noise amplitude tends to oo as the stretched time t' tends to — oo; and 

3. The signal ends at t' — 0. 

As a result, as long as the ratio of the signal amplitude to the noise PSD does not decrease 
too rapidly with decreasing frequency (the precise rate determined by the rate at which 
the frequency changes with time), then p 2 for a inspiralling binary system should increase 
without bound as the observation period extends into the infinite past. This is the role 
that the moment plays in equation |2.4j for p 2 : if the detector noise PSD increases as or 
less rapidly than f~ 7 ^ 3 as / — * 0, then / 7 / 3 diverges as the observation period is extended 
into the infinite past (i.e., as the cut-off frequency /; tends to zero) and the SNR increases 
without bound. In the case of LIGO, the PSD owing to seismic noise increases as f~ A at 
frequencies below the suspension pendulum frequency, so that even an infinite observation 
period leads to a finite SNR. 

Like the SNR, the frequency /o of a truly monochromatic signal (i.e., eqn. |5.20| ) can be 
determined to arbitrary precision given a sufficiently long observation period. By analogy, 
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this is equivalent to the determination of the mass parameter Ai of a inspiralling binary 
system (cf. eqn. |5.23c|) . Consequently, we expect that as long as the ratio of the signal 
amplitude to the noise PSD does not increase too rapidly, the variance in the M. decreases 
to zero as the observation period extends into the infinite past. Too rapidly, in this case, is 
f~ 17 ' 3 . Thus, even though the signal power in a given bandwidth may be much lower than the 
noise power in the same bandwidth, the information present can still play an important role 
in determining the precision with which the parameterization of the signal can be determined. 

Again, we emphasize that these conclusions refer only to the idealized case of a cir- 
cular binary system of two point masses evolving exclusively owing to the emission of 
quadrupole gravitational radiation. The relevance of these conclusions is that the limit 
fi — ■> (i.e., /17/3 — > oo), which may seem far from the reality of observation, is in fact very 
close to that which can be attained in LIGO operating in a regime where all our approxima- 
tions are valid. 

In the limit f)—>-0 (/17/3 — > 00), the variance and the correlation coefficients van- 
ish, corresponding to the determination of M. to infinite precision. The remaining variances 
are (cf. eqns. ggga] , |3T22c1 , |^22d|) 



00 



'I 



P 



-2 



fl 



/3 



fl/3 ~ fi/3 p2 



4/3 



and the remaining correlation coefficients are (cf. eqns. |3.23b| , |3.23c| , |3.23t| ) 



(5.24a) 
(5.24b) 

(5.24c) 



oofrjT — 



f'4/3 

■jl/2- 
J 1/3 



(5.25a) 
(5.25b) 

(5.25c) 



In this limit, the moments fa/3, / 4 / 3 , and f l / 3 describe completely the precision with which 
A, ipi an d T can be measured. For LIGO/VIRGO-like interferometers, these moments 
change negligibly when we pass from fi = 10 Hz to the limit of f\ = Hz [recall that we are 
assuming Sh(f) oc f~ 4 at frequencies below the resonant frequencies of the seismic isolation 
circuit]; consequently, in observing more than the last several minutes of binary inspiral the 
SNR of the observed signal is unchanged, and the variances and correlation coefficients are 
independent of the details of the ultralow frequency behavior of the interferometers. 
By analogy with £^ and S T we define ooS^, and oo^t'- 



0.232 rads initial interferometer, 
0.190 rads advanced interferometer, 

1.32 x 10~ 4 s initial interferometer 
2.71 x 10~ 4 s advanced interferometer, 



(5.26a) 
(5.26b) 
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OO^- J '0 



/ 



1/3 



(/l/3 ^4/3 
/l/3 ~~ ^4/3 




fx 



/3 



/l/3 ^4/3, 



/ c =100 Hz, /; =10 Hz. 



1/2 



/ c =100 Hz,/,=10Hz 



(/l/3 ^4/3 



(5.26c) 



(5.26d) 



The factor oo^r is shown together with in figure ^| for both the initial and advanced 
interferometers. Over the range 50 Hz < f c < 1 Khz the differences between ooSj and £, are 
small. 

Since the SNR is unchanged as the observation period expands from the last several 
minutes of binary inspiral to include the entire lifetime of the binary, observations over more 
than the last several minutes of the lifetime of a binary system will have an insignificant 
effect on the number of binaries observed (cf. eqn. |5.5| ). Comparing the expressions given 
above for with those for <jj (cf. eqns. |3.22a| - |3~22d|) shows that increasing the observation 
period increases the precision with which T can be determined by approximately 10% and 
the precision with which ip can be measured a factor of two. 



VI. CONCLUSIONS 

Inspiralling binary systems of compact objects are regarded as the most certain observ- 
able source of gravitational radiation for the the Laser Interferometer Gravitational-wave 
Observatory (LIGO). As a start toward understanding the capabilities of the LIGO instru- 
ments in observations of inspiralling binary systems, we have investigated the sensitivity 
of a single interferometer of the LIGO type to the gravitational radiation from inspiralling 
binary systems in the quadrupole approximation. 

Observation of binary inspiral in a single LIGO/VIRGO-like interferometer can, in prin- 
ciple, determine a characteristic mass A4, signal amplitude A, time T, and phase ip. The 
mass M. is a function only of the masses of the system's components and its cosmological 
redshift. The amplitude A is inversely proportional to its luminosity distance and depends 
also on a function of four angles describing the relative orientation of the binary and the 
interferometer. The time T is related to the moment of binary coalescence. Finally, if) is 
related to the phase of the binary system at a fixed moment of time and is also a function 
of the relative orientation angles. 

The probability that the detector response is consistent with the presence of a signal from 
an inspiralling binary system is related to the signal-to-noise ratio (SNR) p that characterizes 
the observation. In practice, a threshold po is chosen and we assert that a signal is present 
in the detector output only if p > po- We characterize our uncertainty in the parameters 
ji = {A, Ai, if), T} that describe the detected binary system by defining volumes V(P) in 
parameter space such that ft G V(P) with probability P. 

When p 2 is large, the probability density from which V(P) is constructed is a multivariate 
Gaussian. Consequently, the determination of V(P) is equivalent to the determination of the 
several variance and correlation coefficients that describe the Gaussian. These coefficients 
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in turn describe the statistical uncertainty in the determination of A, Ai, ip, and T, and the 
correlation in the errors in each. 

For observations of binary systems in LIGO/VIRGO-like interferometers, the expected 
SNR, variance, and correlations coefficients may be expressed in terms of the mode of the 
probability distribution P(fJ>) and several moments of the noise PSD of the interferometer 
We have used a detailed model of the PSD for both the initial and advanced LIGO inter- 
ferometers configured for standard recycling, and have evaluated the moments of the PSD, 
the expected SNR, variances, and correlation coefficients as functions of the recycling knee 
frequency 

The two interferometers of the LIGO detector share nearly the same orientation. Con- 
sequently they will act similarly to a single, more sensitive interferometer. In addition to 
providing results for a single interferometer of the LIGO/VIRGO type (either initial or ad- 
vanced), we also express our results for the LIGO two- interferometer network in the limit 
that the interferometers share exactly the same orientation. 

From the expected SNR and an estimate for the cosmological rate density of inspiralling 
binary systems we have calculated the rate of observed binary inspiral events as a function 
of the SNR threshold. We find that for the advanced LIGO detector a conservative estimate 
of the rate of observed binary neutron star inspiral events is 69 yr -1 , of which 7 per year 
will be from binaries at distances greater than 950 Mpc. This is important for observational 
cosmology, since the differential rate (i.e., dN / dAdAi) depends on the Hubble constant and 



other cosmological parameters [16 



For observed binary systems, the fractional standard deviation in the characteristic wave- 
form amplitude is equal to 1/p: if the threshold p is 8, then the fractional one-sigma uncer- 
tainty in the measured amplitude will be less than 12.5% for sources observed in either LIGO 
or LIGO-like interferometers. The chirp mass can be measured to phenomenal precision: 
again, if the threshold p is 8 then the fractional one-sigma uncertainty in At will be less 
than 2.2 x 10~ 5 for binary neutron star systems observed in the advanced LIGO detector. 
We have also calculated the precision with which T and can be determined. 

The optimum detector configuration for the observation of binary inspiral depends sen- 
sitively on the goal of the observation. For example, if the object is to maximize the rate 
of observed binary systems without constraining the uncertainties in A, Ai, T and ip, then 
one detector configuration is clearly favored. On the other hand, if the object is to be 
able to characterize as precisely as possible one of the observables (e.g., Ai), while allow- 
ing that some otherwise observable sources may be missed entirely, then another detector 
configuration is preferred. We have given a concrete formulation to the question of op- 
timum interferometer configuration and answered it in the context of our model for the 
interferometer noise PSD. 

The quadrupole approximation is useful for our LIGO/VIRGO appraisals; however, the 
neglect of post-Newtonian contributions (including spin-orbit and spin-spin interactions) 
to the gravitational radiation luminosity and waveform is a weakness of our estimates and 
should be remedied in a more detailed appraisal of interferometer sensitivity!. Including 



Preliminary Monte Carlo investigations by Cutler |5C] suggest that the inclusion of some of the 
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these interactions will increase the information that can be extracted from gravitational 
radiation observations of binaries over that described here. The formalism that we have 
developed — where the SNR, variances, and correlation coefficients are all expressed in 
terms of moments of the interferometer noise PSD — should prove valuable in that regard: 
it is readily extended to encompass an arbitrarily more sophisticated gravitational radiation 
waveform that is richer in information regarding the source than the one we have studied 
here. 
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terms neglected in our analysis increase the fractional one-sigma uncertainty in M by no more 
than a factor of ten over our estimate, and have a much smaller effect on A, T, and if). 
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FIGURES 



FIG. 1. The noise power spectral density (PSD) Sh(f) for the (anticipated) initial LIGO in- 
terferometers configured for standard recycling with a knee frequency of 300 Hz. The solid line 
shows the total PSD, while the dashed lines show the important physical limits and environmental 
influences that contribute to the total. For more detail, see Sec. |TVA| and table 0. 



FIG. 2. The noise power spectral density (PSD) Sh(f) for the (anticipated) advanced LIGO 
interferometers configured for standard recycling with a knee frequency of 100 Hz. The solid line 
shows the total PSD, while the dashed lines show the important physical limits and environmental 
influences that contribute to the total. In the limit that the two LIGO interferometers have identical 
orientations and we ignore the information available owing to gravitational wave burst arrival time 
differences, the (advanced design) LIGO detector noise PSD Sh(f) is 1/2 the value shown here. 
For more detail, see Sec. |IVA| and table 0. 



FIG. 3. The sensitivity of a LIGO-like interferometer to the gravitational radiation from a 
coalescing binary system depends on the detailed characteristics of the interferometer through 
several moments of the inverse of the its power spectral density (PSD) Sh(f)- In particular, the 
signal-to-noise ratio (SNR) p 2 depends on a moment of S^ 1 (f). Here we show how this moment 
(normalized to its value at a recycling frequency of 100 Hz) varies with the choice of interferometer 
recycling frequency. To maximize the rate at which sources are detected this quantity should be 



maximized. For more details see IV B 



FIG. 4. We define the range function 7£ 7 of a LIGO-like interferometer as the distance di within 
which a fraction 7 of the observable sources are expected to lie. We also define a characteristic 
distance 1Z, such that the total rate of observable sources is A-kTZ^N /3, where M is the rate density of 
sources (which we assume to be uniform). A conservative estimate of 1Z for an advanced LIGO-like 
interferometer is 420 Mpc. Here we show 7 as a function of TZ-y/TZ. For further discussion see 



Sec. VB. 



FIG. 5. The fractional standard deviation a^\jM. of the measured mass M. depends on the 
distance to the source, the relative orientation of the source and the interferometer, and a fac- 
tor YiM(fc) that depends on the interferometer configuration (i.e., the recycling frequency / c ; 
cf. eqn. |5.15| and Sec. VCj ). Here we show £_a/( as a function of f c for initial and advanced 



LIGO-like interferometers. In order to maximize the precision with which A4 can be determined, 
the recycling frequency should be chosen to minimize E^vf. The corresponding recycling frequency 
differs from that which should be chosen to maximize the rate of sources detected (cf. fig. 0). For 



more details, see Sec. V C 
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FIG. 6. The characteristic time T (related to the moment when coalescence occurs) can be 
determined from pre-coalescence observations. The precision with which it can be determined 
depends on (among other things) a factor (defined in Sec. V C| ) that encapsulates all dependence 
of the precision on the duration of the observation and the characteristics of the detector noise. 
The duration of the observation will typically be several minutes, during which time the radiation 
frequency will sweep from 10 Hz to approximately 1 Khz, and assumes that the observation 
period is limited in this way. Here we show Et as a function of the detector recycling frequency for 
both initial and advanced LIGO-like interferometers. An interferometer optimized to measure T 
as precisely as possible should minimize Also shown are factors oo^t for initial and advanced 
interferometers. These are defined identically to except that they assume an observation period 
that extends over the entire life of an idealized binary system evolving only owing to the emission of 
gravitational radiation. The small differences between oo^T and Y,? show that the characteristics 
of the interferometer do not change significantly as the observation period is lengthened beyond 



the last several minutes of binary inspiral. For more details, see Sec. VD 



FIG. 7. The correlation coefficients describe correlations among the statistical errors in the 
measurements of the parameters describing an inspiralling binary system. There is no correlation 
between errors in the measurement of the waveform characteristic amplitude A and any of the other 
measurable parameters. The remaining correlation coefficients are shown here, for both initial and 
advanced LIGO-like interferometers, as a function of the interferometer recycling frequency. For 



more details, see Sec. V C 
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TABLES 



TABLE I. The amplitude of the gravitational radiation waveform observed in a single detector 
depends on the relative orientation between the source and the detector through the function 2 
(cf. eqn. |3.31 ). The orientation angles are unknown and cannot be acertained by observation; 
however, their a priori distribution is known and consequently the a priori distribution of the 
signal amplitude for binaries at a fixed distance from the detector is also known. In this table we 
give the cumulative probability distribution of O 2 , and also a function of P(@ 2 > x 2 ) that arises 
when we evaluate the number of sources within a given distance whose signal amplitudes exceed a 
given threshold. For more details, see Sec. HIE, V A , and [V Bj. 
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TABLE II. The sensitivity of LIGO depends in large measure on the details of the instrumen- 
tation, and this is expected to evolve over the lifetime of the physical plant. We consider the two 
extremes of detector technology that have been cited by the LIGO team in their reports Jl|,0|: 
the instrumentation that is expected to be available when the facilities first come online (initial 
interferometers), and that expected to be available much later (advanced interferometers). 



Initial 


Advanced 


interferometers 


interferometers 



Temperature (T) 


300°K 


300°K 


Pendulum frequency (/o) 


1 Hz 


1 Hz 


Suspension Quality (Qo) 


10 7 


10 9 


End mass (m) 


10 Kg 


1000 Kg 


End mass fundamental mode (/ m t) 


23 KHz 


5 KHz 


End mass quality (Q m t) 


10 5 


10 5 


Effective laser power (lor]) 


5 W 


60 W 


Laser wavelength (A) 


5139A 


5139A 


Mirror losses (A 2 ) 


5 x 10" 5 


2 x 10~ 5 



TABLE III. The signal-to-noise ratio, the variances, and the correlation coefficients all depend 
on the characteristics of the interferometer through the moments of its inverse noise power spectral 



density (c/. Sec. IV A). Here we give those moments for the initial and advanced LIGO instru- 
mentation for the two cases case where the standard recycling "knee" frequency f c is 100 Hz and 
500 Hz. In both cases, the low-frequency cut-off (determined by the duration of the observation) 
is 10 Hz. 
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